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EXPERIMENTAL EVALUATION OF ATMOSPHERIC EFFECTS ON RADIOMETRIC 
MEASUREMENTS USING THE EREP OF SKYLAB 

1 . Program Summary 

At the present time we are awaiting the results of the evaluation of 
our recently submitted proposal to extend the present contract. 

2. Work Accomplished During the Reporting Period 

In addition to the preparation of a proposal in response to a request 
from the Contracting Office (RTAO: BB631 B4), three other tasks were per- 

formed during the reporting period. These include: 

a. Completion of an "interim" final report to review and sum- 
marize the work completed for the year’s period commencing 
7 May 1974. This document was submitted along regular 
channels to the technical monitor. 

b. A comprehensive evaluation of the information received at 
the recently held SKYLAB EREP Principal Investigators Data 
Meeting with particular attention to how this study will 

be effected by problems encountered with the data. Conclu- 
sions drawn from this evaluation are summarized in Section 

3. 

c. Study of the theoretical model adopted to compute atmo- 
spheric effect including a comparison of some preliminary 
calculations with S191 digital data. These results are 
included in Section 4. 

3. Evaluation of Material From the PI Data Meeting 

In a series of meetings and discussions with J. Barnes and J. Willard 
who attended the meeting and through study of the printed material received, 
a comprehensive evaluation of the impact of data/instrumental problems on 
fulfilling the scientific objectives of this study was undertaken. The 
results of this study are summarized by instrument. 



3.1 S190 A Multispectral Camera - S190 B Earth Terrain Camera 


Since data from these photographic instruments is utilized entirely to 
document the areal homogeneity of surface targets, anomolies such as streaks 
electrostatic marking, fogging, and emulsion lifting are of little conse- 
quence. These imperfections do not prevent imagery from being used to seek 
homogenous areas since for the most part, they are readily distinguishable 
against desert-like backgrounds. — 

3.2 S191 Interferometric Spectrometer 

Data from the long wavelength (LWL) portion of the SI 91 Spectrometer 
scan was to have been used to deduce ground-truth data of atmospheric struc- 
ture over the selected target sites. In addition, since the selected test 
site area was generally acquired up to 45° ahead of nadir, an effective 
viewing angle scan was performed by the astronaut as the site was kept within 
the instrument field of view providing a means of examining the angular 
isotropy of surface reflectance (to validate the Lambertian assumption made 
in the theoretical treatment) in the short wavelength (SWL) portion of the 
data. Examination of the instrumental error problems encountered leads us 
to believe that the LWL data may not be suitable for inferring vertical 
atmospheric structure. Absolute accuracies cited are "15% in the SWL and 
thermal regions. 

Of the problems cited in "S 191: Cautionary Note for Data Processed 

According to Photr 524 CH 2 " forwarded to Pi's recently, the most signifi- 
cant to our study is that of off band radiation. Atmospheric vertical 
structure (particularly the water vapor profile) should be obtainable from 
radiance measurements made in the 6-8 ym (6.7 yin H^O Band) and 14-15 ym 
(15 ym CO^ "hot" Bands) . However, such inversion procedures are very sensi- 
tive to noise in the data. Since the offband radiation problem^seems to " 
become critical in precisely these absorption band spectral regions, it does 
not seem feasible to use the data for this purpose. The only recourse is to 
account for site atmospheric parameters using nearby conventional* meteoro- 
logical data. Climatology would probably provide better results than 
extremely noisy inversions. Due to cool down problems, there is no useful 
thermal data for Passes 2 or 5 which include both ground-truth test sites. 



Data in the SWL region is affected by both calibration and off-band 
radiation problems, but appears to be useable in the region from .45 ym 
2.00 yin. Data from an acquired test site (Pass 43, Salton Sea Desert) for 
a period of ~3.7 seconds shows no apparent change in measured radiance with 
viewing angle. (See Figure 4.) However, these plotted data are not correc- 
ted for atmospheric effects. 

3.3 S192 Multi spectral Scanner 

The low frequency and high frequency noise filtering techniques which 
have been developed are certainly more than adequate to satisfy the data 
requirements of this study. In selecting an areal segment of S192 digital 
data for analysis, areal averaging will most likely be done at a resolution 
of “25-35 km over the homogeneous targets selected for study, this averaging 
will provide good characteristic numbers. 

The only serious impact due to S192 anomolies involves the loss of data 
during Pass 2 due to improper spectral alignment and Pass 5 since C/D/P was 
not properly seated. Fortunately, alternate segments (Passes 39 and 43) 
were requested for these sites, but simultaneous ground truth has been lost. 

4. Preliminary Theoretical Results 

Theoretical values for the S192 channel radiances are calculable 
from the expression: 

X 2 

R 1 (a,b, c) = if I(X) ^ W [T 0 W r (X) sin 0 (aT z (X) - b) + cJ Q (X,t)] dX 

Vx\ 

Radiances at zero air mass (such as are measured by aircraft ground ^ 
truth) are similarly calculated by assigning unit transmissivity to the 
atmosphere on the exit path. The difference between the theoretical radiance 
and the zero air mass radiance is the atmospheric correction for the channel, 
where: (Sources and references in brackets) 

1) Superscript i refers to a particular SI 92 channel with wave- 

i i 

length limits specified by (X^, X£) *- i = 1,2,3. . 13. 



2) Parameters (a,b,c) identify the calculated intensity R 1 in 
the following manner: 

R 1 (1,0,1) = theoretical S192 channel radiance 

R 1 (0, -1,0) = Theoretical "zero air mass" radiance 

R 1 (1,1,1) = R 1 (1,0,1) - R 1 (0,-1,0) = atmospheric correction 


3) 1(A) 

4) ^(X) 

5) e 

6) T g (X) 

W 

7) r(X) 

8) J 6 (X,r) 


-2 -1 

solar irradiance spectrum (watts m pm ) 
[Thekaekara, 1970] 
til 

SI 92 i— channel spectral responsivity 
[NASA- JSC, 1973] 

solar elevation angle for the particular test site 
geometry 

atmospheric transmissivity from surface to space at 
angle 6 

zenith atmospheric transmissivity 
[Selby and McClatchey, 1972] 

surface reflectance spectrum 

[ground truth and/or Hunt et al, 1970-1974] 

atmospheric diffuse reflectivity at angle 9 and for 
surface reflectance r 

[Fraser, 1973 and Plass and Kattawar, 1968] 


These computations provide not only for immediate magnitude comparison 
with actual S192 radiances, but allow interchannel (differential) spectral 
effects to be evaluated. 

Simulated S191 data may be obtained by -evaluating only the integrand 
as a function of wavelength while assigning the response function an effec- 
tive value of 1.0 at all wavelengths. 

A sample calculation made with this scheme is illustrated in Figures 1, 

2, 3. Figure 1 lists the input parameters used in the evaluation of the 
equation for radiance. The solar zenith angle is 37.55°. A summer mid- 
latitude atmospheirc model was chosen with a 23 km visual rangle haze model. 
Atmospheric transmissivity (4 § 5) is calculated as a function of wave- 
number (1 52). Columns 3, 6, 7 are input, column 7 being carefully chosen to 



compliment the atmospheric model and viewing geometry. 

Unit band pass (i.e. square response function centered at given wave- 
number) quantities are given in Figure 2. Columns la, lb, lc are, respect- 
ively: theoretically satellite measured radiance, air craft (low flying) 

measured radiance, and atmosphereless radiance. This last quantity represents 
the radiance measured in the total absence of an atmosphere i.e., the 
reflectance spectrum of an identical surface area element illuminated by 
solar radiation in identical viewing geometry at the top of the atmosphere. 
Column 2 is the correction which must be subtracted from a satellite 
measurement to give the radiance which would be measured by an earth 
resources aircraft. Columns 3 b § c are the surface reflectance changes 
corresponding to the atmospheric effect over the zero air mass and atmospheric 
effect over the zero air mass and atmospheric cases, respectively. The larger 
this number is | Abs. Mag. | <_ 1, the greater the atmospheric effect. Some 
of these data (3b) are plotted in Figure 5. Note that the effect of 
Rayleigh scattering is significant, while the aerosols form a continuum. 

Column 4 would be input data of actual measurements if any [here (4)=(la)] 
and Column 5 is the measured data corrected by Column 2 [in this case 
(5) = (lb)]. 

In Figure 3 these results are integrated over the S192 filter 
response functions which the program stores internally. The tabulated 
quantities are analougous to those in Figure 2 differing only in their 
bandpass weighting. The column labelled "GNUBAR" refers to the filter 
function centroid wavenumber as calculated from the stored filter function. 
Here again, the "MEASURED RADIANCE" column is a dummy set equal to the 
theoretical value since no data was available. 

The results of these calculations (Figs. 2 and 3) and plotted in 
Figure 4 and compared to three time segments of S191 data from Pass #3 
over the Salton Sea Desert. Note that the calculated values are too 
large by a factor of 10 (they are plotted as multiplied by 10 *) . The 
reason for this discrepancy is not known. The S191 data was taken from an 
Intermediate Radiance Parameters Tabulation (Product S042-3) . The theoret- 
ical curve is somewhat low in the visible (.4-. 7 ym) , but this is undoubted- 
ly due to the low visible reflection spectrum assumed for the surface (~ .20). 



The column labelled "SURF. REF. CHANGE - ZERO AIR MASS" provides a 
measurement of the atmospheric correction which must be made in surface 
reflectance determinations to interpret space platform - based multi- 
spectral data in times of low altitude aircraft measurements. As expected 
the most significant atmospheric effect is in the visible and slopes down 
into the rear infrared. (Note that since only a few points are plotted, 
the effect of narrow absorption lines is not seen in the data. 

5. Significant Results 

Since the calculations are of a preliminary nature there are no 
significant results to report at this time. 

6. Further Plans 


Further analysis will await receipt of requested data. 
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0 * 1 2 2 3t -02 

0,23586-02 


0,10296-02 

0, 1262E 

00 

-,54996-02 

0,22526-02 

0, 12236-02 

8500, 

1,176 

0.2724E-02 

0, 16176-02 

Q,2Q49t-Q2 


0.11076-02 

0.102/6 

00 

0,49436-01 

0,27246-02 

0,161 /t-02 

8000, 

1,250 

0 * 25 3 3E -02 

G. 14936-02 

0,17666-02 


0.10406-02 

0,10456 

00 

0,65086-01 

0,25536-02 

0 , 1 493t *02 

7500, 

1,313 

0.1320E-02 

0.5629E-03 

0, 1495E-02 


0,73696-03 

0,18966 

00 

-,17506-01 

0,13206-02 

0,58296-03 

7000, 

1,429 

0, 9478E-03 

0,474 7E-04 

0,13016-02 


0 , 90 0 4E-0 3 

0.2645E 

01 

-,40fa8£-01 

0,94/06-03 

0.4/476-04 

6500, 

1,538 

0,165/6-02 

0,9466E-03 

0,10966-02 


0,71036-03 

0.1125E 

00 

0,76856-01 

0,16576-02 

0,94666-05 

6000, 

1,667 

0,13326-02 

0,75326-03 

0,87076-03 


0,57886-03 

0.1155E 

00 

0,79 4 8E-0 1 

0,13326-02 

0,75326-05 

5500, 

1,818 

0,47866.-03 

0*6 1 58E-07 

0,61586-03 


0.4766E-03 

0,11666 

04 

-,33426-01 

0,47866-03 

0,61706-07 

5000, 

2,000 

0.5541E-03 

0,9001 t-04 

0,41546-03 


0 ,26416-03 

0,440l£ 

00 

-.22136-01 

0,35416-03 

0,90016-04 

4500, 

2,222 

0,25176-03 

0, 3093E-Q7 

0.3093E-03 


0,251 7E-03 

0,12206 

04 

2794E-0 l 

0,25176-03 

0,30976-07 

4000, 

2,500 

0 , 1646E-03 

0 ,21746-OS 

0 ,22}86-03 


0, 1&24E-03 

0 , 1 259E 

02 

-,25186-01 

0, 1846E-03 

0 ,2 1 746*05 



3 . 

COMPUTATION UF 3192 CHANNEL WEIGHTED RADIANCES 


CHAN 

NELA 

GNU* 

bar 

CHANfiEL 

LIMITS 

ThEURETICAL 

RADIAN Ct 

ZERO AIR MASS 
RADIANCE 

ATMDSPHERELESS ' 
RADIANCE 

THEURETICAL 

CORRECTION 

SURF ,HEF , CHANGE 
ZERO AIR HASS 

SURF, REF .CHANGE 
NU ATMOSPHERE 

MEASURED 

RADIANCE 

currected 

radiance 

1 

23061, 

20900, 

25000, 

0,13806-01 

0,33/46*02 

0,1 3816*0 1 ' 

0,1043E-01 

0,43576 00 . 

0,1 4556 00 

0,13806-01 

0,33746*02 

2 

21406, 

16600, 

23600, 

0,14276-01 

0,44966-02 

0,14276*01 

0,97696*02 

0,34296 00 

0,1 1 69E 00 

0.14276-01 

0,44966-02 

i 

19592, 

17300, 

21700, 

0.1336E-01 

0 ,53496*02 

0,1 3376-01 

Q,8G13E-02 

0 ,2642E 00 

0,66136-01 

0,1 336E-0 1 

0,53496-02 

0 

16176, 

15900, 

20400, 

0,12196-01 

0 ,56266-02 

0,12196*01 

0,63606*02 

0.2147E 00 

0.6037E-01 

0,12196*0 1 

0.58266*02 

5 

1 6963 | 

15000, 

16500, 

0,11356-01 

0 , 6333t «02 

Q.1136E-01 

0.5020E-02 

0, 1 717E 00 

0,35226*01 

0,11356-0 1 

0,6333E-02 

6 

15268, 

13400, 

16900. 

0,961 QE-02 

0 ,65716*02 

0,96206-02 

0.3239E-0 2 

0.6496E 02 

*,7557E*02 

0 , 98 1 0fc*02 

0.65716-02 

7 

13196, 

10000, 

15100, 

0,71116-02 

0 , 39936*02 

0,7U2E-Q2 

0,31176-02 

0 ,5 1 6 1 E 03 

*, 7651 E-0 1 

0,71116-02 

0.3993E-02 

6 

10739, 

9100, 

12900, 

0,40666-02 

0 , 18586-02 

0,40716*02 

0, 22096*02 

0,40916 03 

- , 3346E-0 1 

.0,40666-02 

0,18586*02 

9 

9726, 

7700, 

11100, 

0.3322E-02 

0 , 19306-02 

0 ,3324E-02 

0,13926-02 

0, 1 718E 02 

0,27676-01 

0,33226-02 

0,19306*02 

10 

6816, 

7200, 

10500, 

0.2721E-02 

0 , 15366*02 

0,2723E*02 

0,1 1656-02 

0 , 1 659E 02 

0,31196-01 

0,27216-02 

0,15366-02 

11 

6396, 

5400. 

7800, 

0.1156E-02 

0,47696-03 

0,11696*02 

0,67666-03 

0,31966 03 

0,23666*0 1 

0,11566-02 

0,47896-03 

12 

4461, 

4000, 

0,35866*03 

0,1420E*03 

0 ,36316*03 

0.2167E-03 

0.3078E 03 

O,2557E-01 

0,35866-03 

0,14206*03 


5200, 
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